Radiating black hole solutions in arbitrary dimensions 
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I. INTRODUCTION 
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We prove a theorem that characterizes a large family of non-static solutions to Einstein equations 
in Ai'-dimensional space-time, representing, in general, spherically symmetric Type II fluid. It is 
shown that the best known Vaidya-based (radiating) black hole solutions to Einstein equations, in 
I both four dimensions (4D) and higher dimensions (HD), are particular cases from this family. The 

, spherically symmetric static black hole solutions for Type I fluid can also be retrieved. A brief 

discussion on the energy conditions, singularities and horizons is provided. 
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' In recent time, it was demonstrated that the string theory requires higher dimensions for its consistency. Models 
with the space-time with large extra dimensions were recently proposed in order to solve the hierarchy problem, that 
is to explain why the gravitational coupling constant is much smaller than the coupling constants of other physical 
I interactions. These new concepts of higher dimensional physics have a number of interesting applications in modern 
^ . cosmology and theory of gravity. This has triggered, during the past decade, a significant increase in interest in black 
p\| ' holes in higher dimensions (HD) (see, for example, the review articles of Horowitz [l[ and Peet [2])- There is now 
, an extensive literature of solutions in string theory with horizons, which represent black holes, and related objects 
in arbitrary dimensions. The physical properties of these solutions have been widely studied. Interest in black holes 
, in HD has been further intensified in recent years, due, for example, to the role they have played in the conjectured 
correspondence between string theory (or supergravity) on asymptotically locally anti-de Sitter backgrounds and the 
] large-N limit of certain conformal field theories defined on the boundary-at-infinity of these backgrounds [1, 0, S] ■ 
QO Static and spherically symmetric space-times are one of the simplest kinds of space-times that one can imagine in 

general relativity. Yet, even in this simple situation, solving the Einstein field equations may be far from trivial. In 
fact, it turns out that such a problem is intractable due to complexity of the Einstein field equations. Hence, there 
are very few inhomogeneous and nonstatic solutions known, one of them is the Vaidya solution. The Vaidya solution 
@ is a solution of Einstein's equations with spherical symmetry for a null fluid (radiation) source (a Type II fiuid) 
' described by energy momentum tensor Tab — i^lah, la being a null vector field. The Vaidya's radiating star metric 
is today commonly used for two purposes: (i) As a testing ground for various formulations of the Cosmic Censorship 
Conjecture (CCC). (Actually CCC is a famous conjecture, first formulated by Penrose The conjecture, in it's 
weak version, essentially state that any naked singularity which is created by evolution of regular initial data will be 
shielded from the external view by an event horizon. According to the strong version of the CCC, naked singularities 
are never produced, which in the precise mathematical terms demands that space-time should be globally hyperbolic.) 
(ii) As an exterior solution for models of objects consisting of heat-conducting matter. Recently, it has also proved to 
be useful in the study of Hawking radiation, the process of black- hole evaporation Q , and in the stochastic gravity 
program Q . It has also advantage of allowing a study of the dynamical evolution of horizon associated with a radiating 
black hole. 

Also, several solutions in which the source is a mixture of a perfect fluid and null radiation have been obtained in 
later years This includes the Bonnor- Vaidya solution [ll| for the charge case, the Husain solution with an 
equation of state P = kp. Glass and Krisch [l^ further generalized the Vaidya solution to include a string fluid, while 
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charged strange quark fluid (SQM) together with the Vaidya null radiation has been obtained by Harko and Cheng 
(see also p^). Wang and Wu [IGj further extrapolated the Vaidya solution to more general case, which include 
a large family of known solutions. 

Motivated by this and by a recent work Salgado fl3| , we [l^ have proved a theorem characterizing a three parameter 
family of solutions, representing, in general, spherically symmetric Type II fluid that includes most of the known 
solutions to Einstein field equations. This is done by imposing certain conditions on the energy momentum tensor 
(EMT) (see also [H [H HH, 113). 

In this paper, we consider an extension of our work [18], so that a large family of exact spherically symmetric Type 
II fluid solutions, in arbitrary dimensions, are possible, including it's generalization to asymptotically de Sitter/anti-de 
Sitter. 



II. THE RADIATING BLACK-HOLE SOLUTIONS 



Theorem - I; Let (M, gab) be an N-dimensional space-time [sign{gab) = +{N — 2)] such that (i) It is non-static and 
spherically symmetric, (ii) it satisfies Einstein field equations, (Hi) in the Eddington-Bondi coordinates where 



ds' 



~A{v,ryf{v,r) dv^ + 2eA{v,r) dv dr + r^{dnN^2) , where, (df^Ar_2)" 



sin^{ei)sin^{e2)del- 



IN-2 



ni- - 
■ -, sm" 



del 
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, the energy-momentum tensor T"^^ satisfies the conditions — 

and Tg^ = kT^^ , {k = const. G TZ) (iv) it possesses a regular Killing horizon or a regular origin. Then the metric of 
the .space-time is given by 



ds' 



1 



2m{v, r) 
{N - 3)r^-3 



dv^ + 2edvdr + r'^{dilN-2)'^ , 



(6 = ±1) 
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M{v) - 87rC(w) ( 
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N-2 J {N-2)k+l ' 
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if C'iv) = 0, 

if C{v) ^ and fc 7^ -l/(^ - 2), 
if C{v) 7^ and fc = -1/{N - 2). 
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^if-2)U^ '^^^S[lA.k,...,k]. (3) 

and 

1^ f N^) dM 1 ^^(jV-2)(fc^l) + l if ^ _ 2), 

Sttv^^ ^ \N — 3 J dv [N --2)k-\-l ov / / \ / ^ 

K - { (4) 
d-- (#51) ^ - Inr if fc = -1/{N 2). 

Here, M(v) and C(v) are the arbitrary functions whose values depend on the boundary conditions and the fundamental 
constants of the underlying matter. 

Proof: Expressed in terms of Eddington coordinate, the metric of general spherically symmetric space-time in A^- 

dimensional space-times (n, H m '261 IS, 

ds^ ^ ~A{v,rff{v,r) dv"^ + 2eA{v,r) dv dr -\- r'^{dflN-2)'^ . (5) 

Here A(v,r) is an arbitrary function. It is useful to introduce a local mass function m{v,r) defined by f{v,r) = 
1 — 2m{v,r)/{N — 3)r'-^~^\ For m{v,r) — m{v) and A = 1, the metric reduces to the A^-dimensional Vaidya metric 

El. 

In the static limit, this metric can be obtained from the metric in the usual, spherically symmetric form, 

ds^ = -f{r) dt^ + — + r\dnN-2?. (6) 
/(n 
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by the coordinate transformation 



dv = Air)-^ {dt + 

fir) 



(7) 



In case of spherical symmetry, even when f{r) is replaced by f{t^r)^ one can cast the metric in the form ([5]) [27 
The non- vanishing components of the Einstein tensor [23 | are 
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where = {u, r, 6'i,... 0^-2}- We shall consider the special case = (hypothesis), which means from 

Eq. ([5a|l . A{v,r) = g{v). This also implies that = G^. (Eq. (j8dp V However, by introducing another null coordinate 
V = J g{v)dv, we can always set without the loss of generality, A{v,r) = 1. Hence, the metric takes the form, 



ds^ = - 



1 



2m{v, r) 
{N - 3)r(^-3) 



dv^ + 2edvdr + r^{dClN-2f 



(9) 



Therefore the entire family of solutions we are searching for is determined by a single function m{v,r). Henceforth, 
we adopt here a method similar to Salgado 17[ which we modify here to accommodate the non static case. In what 
follows, we shall consider e = 1. The Einstein field equations are 



(10) 



and combining Eqs. ([8]) and (|10p . we have if a ^ 6, T§ = except for a non-zero off-diagonal components T^. In 
addition, we observe that — TJ. Thus the EMT can be written as : 



Tia 
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which in general belongs to a Type II fluid with T^l ^Tg'- = ... = Tg^Zl I* may be recalled that EMT of a Type II 
fluid has a double null eigen vector, whereas an EMT of a Type I fluid has only one time- like eigen vector [2^. On 
the other hand, from the Einstein equations, it follows that 



Enforcing the conservation laws VaTj^ ~ 0, yields the following non-trivial differential equations: 
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and, using, = , 
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Using the assumption made above that Tg^ = kT^, we obtain the following linear differential equation 

5t; (N - 2) 



dr 



which can be easily integrated to give 



C{v) 



j,(iV-2)(l-fe) ' 



(14) 
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where C{v) is an arbitrary function of v, arising as an integration constant. Then, using hypothesis (iii), we conclude 
that 



Civ 



YrfeTdiag[l,l,fc, 



Now using Eqs. ^ [with f{v,r) = 1 - 2m{v,r)/{N - ^ and we get 
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Civ) 



N-2j 7.-(^-2)fc ' 
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which trivially integrates to 
( M{v) 

Af(z;)~87rC(«)(f5f 



m(w, r) 
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^ M{v)~%nC{v)[^) 



(N-2)k+l 
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if Civ) = 0, 

if Civ) 7^ and fc 7^ -l/(^ - 2), 
if Civ) 7^ and fc = -1/iN - 2). 



(18) 



Here the function Miv) arises as a result of integration. What remains to be calculated is the only non-zero off-diagonal 
component of the EMT. From Eqs. ^ and HU]), one gets 



r|-^r 



1 



iV-3 



dm 
dv ' 



(19) 



which, on using Eq. ([T51). gives 

J ( N-2 \ dM 

r"-2 yN-3 J dv 



{N-2)k+l dv ' 



ac (N-2)(k-l) + l 



1 ( N-2 \ dM 1 ac i„ , 



iffc^-l/(7V-2), 
if fc = -l/(7V--2). 



(20) 



It is seen that Eq. is identically satisfied. Hence the theorem is proved. The theorem proved above represents a 
general class of non-static, A^-dimensional spherically symmetric solutions to Einstein's equations describing radiating 
black-holes with the EMT which satisfies the conditions in accordance with hypothesis (iii). The solutions generated 
here highly rely on the assumption (im). On the other hand, although hypothesis (iv) is not used a priori for proving 
the result, but it is indeed suggested by regularity of the solution at the origin, from which, — T^\r=o (see [l3| for 
further details). 

The family of the iV-dimensional solutions outlined here contains A^-dimensional version of, for instance, Vaidya 
[H, |29^ Bonnor- Vaidya [H, H^, dS/AdS [2^, global monopole [M HI, [Ml, Husain [H, HI El, and Harko-Cheng 
SQM solution 14, 23]. Obviously, by proper choice of the functions M(i>) and Civ), and fc— index, one can generate 
as many solutions as required. The above solutions include most of the known Vadya-based spherically symmetric 
solutions of the Einstein field equations, ^^hen N = 4, the 4Z) solutions derived in [l^, [l^ can be recovered. The 
static black holes solutions, in both HD [3_2] and in 4D [l3|, can be recovered by setting Miv) = M, Civ) = C , with 
M and C as constants, in which case matter is Type I. 

In summary, we have shown that the metric 



2Miv) 



WnCiv) 



^(Af-2)(fc-l)+2 



(7V-3)r^-3 iN - 2)[iN ~ 2)k + I] 



dv"^ + 2dvdr + r'^idVLN^2)'^ 
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is a solution of the Einstein equations for the stress energy tensor Eqs. ([3]) and (|4]). A metric is considered to be 
asymptoticahy flat if in the vicinity of a spaceUke hypersurface its components behave as 

, <^ab{xVr,t) , 1 \ 

gab riab H + O -r-- , (22) 



as r — > oo. (e > 0, rjat is the Minkowski metric, aat is an arbitrary symmetric tensor, and is a flat coordinate system 
at spacehke inflnity). According to this deflnition, our metrics Eq. ([1]) are asymptotically flat for k > —1/{N — 2) and 
are cosmological for k < —1/{N — 2). In particular, for k=-l, M{v) = M and 2C{v) = Q^, the metric is just higher 
dimensional Reissner-Nordstrom. The detailed of the asymptotic structure of spatial infinity in higher-dimensional 
space-times can be found in Ref. [ssj and the different conformal diagrams for maximal extension of 4D Vaidya is 
discussed by Fayos et al. [33 |. 

The theorem proved shows that rather than a mathematical coincidence, the above form of the metric is a conse- 
quence of the features of the energy-momentum tensor considered. In the above exact solutions, the associated energy 
momentum tensors share some properties that are taken into account in the theorem in a general fashion without 
specifying the nature of the matter. Therefore, the theorem helps to characterize a whole two-parameter family of 
solutions to the Einstein field equations. 

The solutions discussed in the section are characterized by two arbitrary functions M(v) and C{v), and the 
cosmological constant A. Thus one would like to generalize the above theorem to include A. We can show that 
the energy momentum tensor components, in general, can be written as, T^" = T^f^i, — -^^t l2^ . where A is 
the cosmological constant and T^f^^^ is energy momentum tensor of the matter fields that satisfy T^f-^g = kTJ^^y. 
A trivial extension of the theorem allows one to cover a three-parameter family of solutions, with one of the pa- 
rameters being a cosmological constant A. Next, we just state (proof being similar) the generalization of the Theorem I. 

Theorem - II.' Let (M, gab) be an N- dimensional space-time [sign{gab) — +{N — 2)/ such that (i) It is non-static 
and spherically symmetric, (ii) it satisfies Einstein field equations, (Hi) the total energy-momentum tensor is given 
by = T^f-^i, ~ ^'^6 ; where A is the cosmological constant and T^^^^^ is energy momentum tensor of the matter fields, 
(iv) in the Eddington coordinates where ds"^ — — A{v , r)"^ f {v , r) dv"^ -\- 2A{v,r) dv dr -\- r^(df2jv-2)^, the EMT T^^f^j, 
satisfies the conditions T^^jy — 0, '^[f)o^ — ^"^{jy ~ const. G TZ), (v) it possesses a regular Killing horizon or a 
regular origin. Then the metric of the space-time is given by metric (QP, where 

M{v) + (iv^Sp^Ar^-i if C{v) = 0, 



m{v, r) = < 



Miv) - SnCiv) (f5f ) + ^^Ml_^Ar^-^ if dv) ^ and M -l/iN - 2), 

Miv) - 8t:C{v) (fei) Inr + (jv-^K^-i) Ar^~' if C{v) ^ and k = -l/(iV - 2). 

(23) 

= ^rp^Stfe) diag[l, 1, fc, . . . , /c] - ^diag[l, 1,1,...,!] (24) 

d-- (f5f ) ^ - W^.^ri--^'^'-''^^ if k ^ -1/iN - 2), 

T: = { (25) 
d-- (^) ^ - if - 2). 

Here, M(v) and C(v) are arbitrary functions of v, arising as integration constants, whose values depend on the 
boundary conditions and the fundamental constants of the underlying matter. 



and 



III. ENERGY CONDITIONS 



The family of solutions discussed here, in general, belon gs t o Type II fiuid defined in [23|. When m = m{r), we 
have /x=0, and the matter field degenerates to type I fluid [16!]. In the rest frame associated with the observer, the 
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energy-density of the matter will be given by (assuming A = 0), 

M = r:, P--r/^-T;^- ^^/g_,^ , (26) 

and the principal pressures are Pi = T? (no sum convention). Therefore Pr ~ = —p and Pgi = kP^ = —kp 
(hypothesis {in)). 

a) The weak energy conditions (WEC): The energy momentum tensor obeys inequality Tabw'^w^ > for any timelike 
vector, i.e., 

/i > 0, p > 0, Pe,=Pe,=... = > 0. (27) 

We say that strong energy condition (SEC), holds for Type II fluid if, Eq. ([27ll is true., i.e., both WEC and SEC, for 
a Type II fluid, are identical. 

b) The dominant energy conditions : For any timelike vector Wa, T°^WaWh > 0, and T'^^Wa is non-spacelike vector, 
i.e., 

/i>0, p>Pei,Pe. =... =/'9,„_„ >0. (28) 

Clearly, (a) is satisfied if C{v) < 0,k < 0. However, p, > gives the restriction on the choice of the functions M{v) 
and C{v). From Eq. (|3]), (fc ^ — l/(n — 2)), we observe p > requires, 

1 [^] ^ _ 1 dC 

8TTrN~'^\N-3j dv {N-2)k + ldv ^ ' 

This, in general, is satisfied, if 

N -2\ dM and, either ^ > and fc < ~1/{N - 2), or ^ < and k > -1/{N - 2). (30) 



N — 3 J dv dv dv 



On the other hand, for k = -1/{N - 2), ^ > if dM/dC > 8n (fff) Inr. The DEC holds if C{v) < and 
— 1 < fc < 0, and the function M is subject to the condition (|30p . Clearly, < — fc < 1. 



IV. SINGULARITY AND HORIZONS 



The invariants are regular everywhere except at the origin r = 0, where they diverge. Hence, the space-time has 
the scalar polynomial singularity [2^ at r = 0. The nature (a naked singularity or a black hole) of the singularity 
can be characterized by the existence of radial null geodesies emerging from the singularity. The singularity is at 
least locally naked if there exist such geodesies, and if no such geodesies exist, it is a black hole. The study of causal 
structure of the space-time is beyond the scope of this paper and will be discussed elsewhere. 

In order to further discuss the physical nature of our solutions, we introduce their kinematical parameters. Following 
York [35] a null-vector decomposition of the metric H]) is made of the form 

gab = -nJb - Lnt + Jab, (31) 

where, 

lab = rX'St^+r' 

= nan'' = /a"" = -1, 

lab = 0; Jab = 0, (32c) 



2m{v, r) 
(N - 3)r^-^ 



(32a) 



Oa K , 



(32b) 
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with m{v, r) given by Eq. ([2]). The optical behavior of null geodesies congruences is governed by the Carter [36!| form 
of the Raychaudhuri equation 



dv 



= /Ce - Rabn" - (7c)" - <^ab<J 



ah I ah 
UJah^^ , 



(33) 



with expansion Q, twist lo, shear cr, and surface gravity IC. Here Rah is the A^-dimensional Ricci tensor, 7^ is the trace 
of the projection tensor for null geodesies. The expansion of the null rays parameterized by v is given by 



e^Va^-JC, 

where the V is the covariant derivative. In the present case, a = ui — [35|], and the surface gravity is, 

IC = -n'^l^Vbla- 



(34) 



(35) 



Spherically symmetric irrotational space-times, such as under consideration, are vorticity and shear free. The structure 
and dynamics of the horizons are then only dependent on the expansion, Q. As demonstrated by York 35], horizons 
can be obtained by noting that (i) apparent horizons are defined as surface such that c± and (ii) event horizons 
are surfaces such that d<d/dv ~ 0. Substituting Eqs. ([23]), ([32]) and into Eq. ([Ml), we get, [k ^ -1/2) 



where, 



e 



T7/ ^ ^(«) 

M{v) - ^ ' 



2M(v 



rN-3 



2 2 

X r 



2A 



(A^-3) 



,X 



{N -2){N -1) 



167rC(w) 



(iV-2)((iV-2)fc+l)' 



(36) 



(37) 



Since the York conditions require that at apparent horizons 8 vanish, it follows form the Eq. ([M]) that apparent 
horizons will satisfy 



(38) 



which in general has two positive solutions. For X'^ — 0^ = 0; we have Schwarzschild horizon r — (2M) "-^ , and 
for M = = we have de Sitter horizon r = 1/x- As mentioned above, for k = —1, one gets Bonnor-Vaidya 
solution, in which case the various horizons are identified and analyzed by Mallett [37l | and hence, to conserve space, 
we shall avoid the repetition of same. For general fc, as it stands, Eq. (|38p will not admit simple closed form solutions. 
However, for 



{{N~2)k + 1) 



2M{{N ~2)k + l) 



with = 0, the two roots of the Eq. 



{N -3){{N -2)k-{N 
coincide and there is only one horizon 

2M{{N -2)k + l) 



(W-4)-(W-2)fc 
(N-3) 



(iV-3)((iV-2)fc- (7V-4)) 



(39) 



(40) 



For < there are two horizons, namely a cosmological horizon and a black hole horizon. On the other hand if, 
the inequality is reversed, > no horizon would form. 



V. CONCLUDING REMARKS 



In the study of the Einstein equations in the 4D space-time several powerful mathematical tools were developed, 
based on the space-time symmetry, algebraical structure of space-time, internal symmetry and solution generation 
technique, global analysis, and so on. It would be interesting how to develop some of these methods to higher 
dimensional space-time. With this as motivation, plus the fact that exact solutions are always desirable and valuable, 
we have extended to higher dimensional space-time, a recent theorem fl8| and it's trivial extension (that includes 
cosmological term A), which, with certain restrictions on the EMT, characterizes a large family of radiating black hole 
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solutions in N-dimensions, representing, in general, spherically symmetric Type II fluid. In particular, the monopole- 
de Sitter-charged Vaidya and Husain solutions can be generated form our analysis and when n = 4, one recovers the 
4D black solutions. If M = C = constant, we have /i=0, and the matter field degenerates to type I fluid and we can 
generate static black hole solutions by a proper choice of these constants. Since many known solutions are identified 
as particular case of this family and hence it would be interesting to ask whether there exist realistic matter that 
follows the restrictions of the theorem that would generate a new black hole solution. 

The solutions depend on one parameter fc, and two arbitrary functions M{v) and C(v) (modulo energy conditions). 
It is possible to generate various solutions by proper choice of these functions and parameter k. Further, most of the 
known static spherically symmetric black hole solutions, in 4D and HD, can be recovered from our analysis. 

It should be interesting to apply these metrics to study the gravitation collapse and naked singularities formation. 
Finally, the result obtained would also be relevant in the context of string theory which is often said to be next "theory 
of everything" and in the study of gravitational collapse. 



Acknowledgments 

Authors are grateful to the referee for constructive criticism. One of the author(SGG) would like to thank lUCAA, 
Pune for hospitality while this work was done. 



[1] Horowitz, G. T. (1992). hep-th/9210119 The dark side of string theory: black holes and black strings (Trieste Spring School 

on String Theory and Quantum Gravity). 
[2] Peet, A. (1999). |hep-th/ 0008241 Lectures on Black Holes in String Theory (TASI). 
[3] Aharony, O., Gubser, S., Maldacena, J., Ooguri, H. and Oz, Y. (2000). Phys. Rep. 323, 184. 
[4] Maldacena, J. (1998). Adv. Theor. Math. Phys. 2, 231. 
[5] Witten, E. (1998). Adv. Theor. Math. Phys. 2, 253. 

[6] Vaidya, P. C. (1951). Proc. Indian Acad. Set. A33, 264; (1999). Gen. Rel. Grav. 31, 119. 

[7] Penrose, R. (1969). Riv Nuovo Gimento 1, 252 ; General Relativity, an Emstem Gentenary Volume, edited by S. W. Hawking 

and W. Israel (Cambridge University Press, Cambridge). 
[8] Parentini, R. (2001). Phys. Rev. D 63, 041503. 

[9] Hu, B. L. and Verdaguer, E. (2004). Living Rev. Relativity 7, 3 htt p : / / www .livingreviews.org/lrr-2004-3'4. 
[10] Krasianski , A. (1997). Inhomogeneous Cosmological Models (Cambridge University Press, Cambridge). 
[11] Bonnor, W. B. and Vaidya, P. C. (1970). Gen. Rel. Grav. 1, 159. 
[12] Husain, V. (1996). Phys. Rev. D 53, 1759. 

[13] Glass, E. N. and Krisch, J. P. (1998). Phys. Rev. D 57, 5945 ; (1999). Glass. Quantum Grav. 16, 1175; (1999). J. Math. 
Phys. 40, 4056. 

[14] Harko, T. and Chang, K. S. (2000). Phys. Lett. A 266, 249. 

[15] Ghosh, S. G. and Dadhich, N. (2002). Phys. Rev. D 64, 127502; Dadhich, N. and Ghosh, S. G. (2003). Gen. Rel. Grav. 
35, 359. 

[16] Wang, A. and Yu, Y. (1999). Gen. Rel. Grav. 31, 107. 

[17] Salgado, M. (2003). Class. Quantum Grav. 20, 4551. 

[18] Dawood, A. K. and Ghosh, S. G. (2003). Phys. Rev. D 70, 104010. 

[19] Kiselev, V. V. (2003). Glass. Quantum Grav. 20, 1187. 

[20] Giambo, R. (2002). Glass. Quantum Grav. 19, 4399. 

[21] Dymnikova, I. (2002). Class. Quantum Grav. 19, 725. 

[22] Gallo, E. (2004). Gen. Rel. Grav. 36, 1463. 

[23] Ghosh, S. G. and Dadhich, N. (2001). Phys. Rev. D 64, 047501. 

[24] Patel, L. K. and Dadhich, N. gr-qc/9909068 

[25] Rocha, J. F. V. (2002). Int. J. Mod. Phys. D 11, 113. 

[26] Barrabes, C. and Israel, W. (1991). Phys. Rev. D 43, 1129. 

[27] Nielsen, A. B. and Visser, M. (2006). Glass. Quantum Grav. 23, 4637. 

[28] Hawking, S. W. and Ellis, G. F. R. (1973). The Large Scale Structure of Space-time (Cambridge University Press, Cam- 
bridge) . 

[29] Iyer, B. R. and Vishveshwara, C. V. (1989). Pramana-J. Phys. 32, 749. 
[30] Chatterjee, S., Bhui, B. and Banerjee, A. (1990). J. Math. Phys. 31, 2208. 
[31] Barriola, M. and Vilenkin, A. (1989). Phys. Rev. Lett. 63, 341. 

[32] Myers, R. C. and Perry, M. J. (1986). Annals of Phys. 172, 304; (1988). Dianyan, Xu. Glass. Quantum Grav. 5, 871; ( 

2002). Das, A. and DeBenedictis, A. (2002). Progress of Theo. Phy. 108, 119. 
[33] Shiromizu, T. and Tomizawa, S. (2004). Phys. Rev. D 69, 104012. 



9 



[34] Fayos, F., Martin-Prats, M. and Scnovilla, J. M. M. (1995). Class. Quantum Grav. 12, 2565. 

[35] York, J. W. Jr. (1984). Quantum Theory of Gravity: Essays in Honor of Sixtieth Birthday of Bryce S. DeWitt, edited by 
S. Christcnsen (Hilgcr, Bristol). 

[36] B. Carter, in General Relativity, edited by S. W. Hawking and I. Israel (Cambridge University Press, Cambridge, England, 
1979). 

[37] Mallett, R. L. (1986). Phys. Rev. D 33, 2201; B.D. Koberlein and R.L. Mallett, Phys. Rev. D 49, 5111 (1994). 



